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5.  Introduction 

At  the  present  time  the  best  way  to  fight  breast  cancer  is  through  early  detection.  The 
size  of  malignancies  in  the  early  stage  is  very  small,  and  the  absorption  of  x  rays  by  the 
malignancies  and  the  breast  tissue  is  similar.  Thus  high  contrast  and  good  resolution  are 
especially  important  in  mammography.  In  current  digital  mammography  systems,  contrast  loss 
occurs  from  scattered  x  rays  that  strike  the  detector  from  all  directions,  and  resolution  is  limited 
by  the  pixel  size  of  the  digital  detector  for  computed  radiography  phosphor  plates  (CR  plates), 
or  the  conversion  efficiency  and  total  number  of  pixels  for  visible  light  CCD  systems. 

Capillary  x-ray  optics  use  total  external  reflection  to  guide  x  rays  through  tiny  hollow 
glass  tubes,  similar  to  the  way  conventional  fiber  optics  transmit  light.  A  capillary  lens  is  a 
bundle  of  hollow  glass  tubes  with  channel  sizes  as  small  as  few  microns.  The  critical  angle  for 
total  external  reflection  is  as  small  as  a  few  mrad  (1.5mrad  for  20  keV  photons).1,2  A  capillary 
lens  can  provide  almost  total  scatter  rejection.  A  well-designed  lens  can  also  provide 
magnification  without  image  degradation  from  a  finite  source  spot  size.  Thus  the  effective 
resolution  can  be  improved  while  a  CR  plate  is  used.  Capillary  optics  can  also  demagnify  the 
image  to  match  a  direct  area  detector,  such  as  a  CZT  detector^4,  which  maybe  available  in  the 

near  future.  Unlike  the  fiber-optic-CCD  combination,  this 
system  does  not  involve  the  use  of  a  phosphor  to 
convert  the  x  rays  to  visual  light.  The  solid  state  detector 
can  provide  almost  1 00%  quantum  efficiency. 


Capillary  optical  mammography  has  been 
cooperatively  investigated  by  Center  for  X-ray  Optics  at 
University  at  Albany  and  the  Radiology  Dept,  of 
University  of  Wisconsin  for  many  years.  Previous  work 
has  demonstrated  the  bright  future  of  capillary  optics  in 
digital  mammography  with  both  contrast  and  resolution 
improvement.5,6  Among  the  problems  found  in  the 
previous  work  are  the  small  size  of  the  lens  and  non- 
uniform  transmission.  Artifacts  caused  by  non-uniformity 
of  transmission  were  observed  in  the  image.  Eliminating 
artifacts  by  image  processing  may  be  a  better  and  more 
economic  way  than  fabricating  an  ideal  lens.  Low 
reproducibility  can  be  a  problem  in  lens  manufacture.  It  is  extremely  important  to  design  a  large 
lens  by  combining  several  small  pieces.  So  a  better  understanding  of  lens  behavior  is 
necessary.  Theoretical  simulation  has  been  proved  to  be  very  helpful  in  understanding  the 
capillary  and  lens  behavior,  and  provide  feedback  to  the  manufacturing  process. Error!  Bookmarknot 
defined.  ^  yyouid  a|S0  pe  easier  to  manufacture  large  lenses  if  they  were  shorter;  for  that  reason, 
lead  glass  optics,  which  would  stop  scattered  radiation  with  a  shorter  length,  were  investigated. 

The  proposed  work  is  focused  on  lens  meaurements,  computer  simulations  and  image 
processing.  During  the  first  year  of  the  work,  a  lens  simulation  has  been  developed  and  applied  to  the 
experimental  data.  The  simulation  has  proved  to  be  helpful  to  understand  the  lens  behavior.  MTF  is 
usually  measured  by  taking  an  image  of  a  slit.  It  is  difficult  to  measure  a  stationary  lens  with  this  method 
when  the  transmission  of  the  lens  is  not  uniform.  A  method  to  measure  MTF  by  taking  an  edge  image 
was  developed  with  background  deduction  technique.  The  lens  simulation  was  expanded  in  the  second 
year,  and  lead  glass  optics  were  investigated.  During  the  second  year,  the  principal  investigator  was 
changed  from  Dr.  Lei  Wang  to  Dr.  Hui  Wang. 


Figure  1  a)  Image  bluring  in 
normal  airgap  magnification,  b) 
Capillary  lens  between  the  patient 
and  detector  eliminates  the  focus 
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Figure  2.  Experimental  setup. 


Figure  3.  Electronic  system. 


6.  Body 

6.1.  DEFECT  SIMULATION 

An  important  step  in  optics  development  was  the  breakthrough  in  modeling 
capacity,  which  has  led  to  an  unprecedented  level  of  theoretical  understanding  of  the 
basic  properties  of  polycapillary  fibers,  especially  fiber  defects.  To  evaluate  the 
experimental  performance  of  polycapillary  fibers,  and  design  capillary  optics,  it  is 
necessary  to  be  able  to  predict  theoretical  behavior  for  complex  geometries.  The 
modeling  program  for  single  fibers  is  based  on  a  Monte  Carlo  simulation  of  simple 
geometrical  optics.  The  computational  speed  is  greatly  enhanced  by  a  reduction  to  two 
dimensions  by  projecting  the  trajectory  onto  the  local  fiber  cross-section.7  Reflectivities 
are  computed  from  standard  tables.8  Significant  recent  progress  has  been  made  in 
understanding  the  effect  of  capillary  profile  error,  waviness,  and  roughness  on  the 
transmission  spectra.9  This  is  extremely  important  in  providing  feedback  to  the 
manufacturing  process.  It  has  been  found  that  extremely  good  fits  can  be  produced  with 
only  two  fitting  parameters  if  a  more  physical  model  of  waviness  is  employed.10 


The  experimental  arrangement  for  single  capillary  measurements  is  shown  in 
Figure  2  and  Figure  4.  An  optical  rail  affixed  to  an  optical  table  carries  an  x-ray  source, 
fiber  platform,  and  x-ray  detector.  Each  can  be  positioned  independently  in  three 


Type 

Description 

Outer  Diameter, 
mm 

Channel  Size, 
pm 

Open  Area  Length, 
mm 

■ 

Borosilicate 

0.5 

12 

65% 

105 

B 

Lead  glass 

0.5 

12 

52% 

95 

Borosilicate 

0.75 

22 

50% 

136 

D 

Borosilicate 

4 

12 

55% 

130 

E 

Borosilicate 

0.3 

4-5 

55% 

105 

Table  1.  Description  of  polycapillary  fibers. 
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■  exp,  Fiber  A  . sim,  Fiber  A 

A  exp,  Fiber  B 

▼  exp,  Fiber  C  . sim,  Fiber 

+  exp,  Fiber  D  - sim,  Fiber 

*  exp,  Fiber  E  - sim,  Fiber  E 


dimensions.  A  collimator  is  placed 
before  the  fiber  and  any  remaining 
x-ray  leakage  around  the  fiber  is 
eliminated  with  metal  powder  or 
filings.  The  measured  fibers  are 
described  in  Table  1. 


30 '  . •  >*<•...  The  results  of  transmission 

y'  studies  as  a  function  of  photon  energy 

20 '  ,  are  shown  in  Figure  4. 11,12  All  of  the 

10.  ■  jjj""*  fibers  except  the  lead  glass  have 

transmissions  at  20  keV  nearly  equal 

0  ^  T  ^  t  ^  ^  to  their  fractional  open  area  (the 

Photon  energy  (KeV)  fraction  of  the  cross  section  of  the 

,  ,  capillary  which  is  open  space,  the  rest 

Figure  4.  Measured  transmission  versus  energy  for  ^  wg  ^  transmission 

polycapillary  fibers  listed  in  Table  1.  .  .  ..  . 

^  y  corresponds  to  the  primary 

transmission  expected  for  a  linear  capillary  optic  employed  as  an  antiscatter  grid. 

The  transmission  falls  off  at  higher  energies.  Modeling  this  transmission  spectra  has 
lead  to  a  greater  understanding  of  defects  in  polycapillary  optics. 

In  Figure  5,  simulations  with  or  without  roughness  corrections  are  compared  with 
the  experimental  data.  In  these  measurements  the  source  is  scanned  transverse  to  the 
fiber  axis.  The  simulation  with  a  roughness  height  of  0.5  nm  fits  the  experimental  data 
quite  well.  It  is  definitely  over-corrected  when  the  roughness  is  1 .0  nm.  The  same 
simulations,  shown  in  Figure  7,  are  also  carried  out  at  68  keV,  where  the  width  of  the 
curve  is  narrower  than  that  at  10  keV  because  of  the  smaller  critical  angle.  As  a  result, 


Source  position(mm) 

Figure  5.  Source  scans  at  10 
keV. 


Figure  6.  Transmission  spectra  of  Fiber  3  simulated 
with  different  bending  curvature  alone  and  compared 
with  experimental  data 
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■  data 

_ z=0 

_ z=1  nm 

_ z=2  nm 


Source  position(mm) 

Figure  7.  Source  scans 
at  68  ke  V. 


photons  also 
experience  fewer 
bounces  on 
average. 
Simulations  with 
roughnesses  as 
large  as  1 .0  and 
2.0  nm  still  could 
not  fit  the  data. 
Knowing  that  1 .0 
or  2.0  nm 
roughness  is 
definitely  too 
large  at  1 0  keV, 
we  can 

determine  that 
the  roughness 
correction  by 
itself  is  not 


Fhctonenergy(KfeV) 


Figure  8.  Simulations  of  transmission  spectra  with 
waviness  only  compared  with  the  experimental  data. 


sufficient  at  high  energy  to  reproduce  source  scan 
measurements.  Other  effects  need  to  be  considered.  These  are  bending  and 


waviness. 


Source  position(mm)  Source  position(mm)  Source  position(mm)  Source  position(mn) 

Figure  9.  Simulated  source  scan  curves  compared  with  experimental  data  at  four  different  photon  energies. 
Parameters  are:  R  —  125  m,  Admax—  0.35  mrad,  roughness  height=0.5  nm. 
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A  slight  bending  of  the 
capillary  can  dramatically  reduce 
the  transmission  of  high  energy 
photons  because  of  the  small 
critical  angle.  A  comparison 
between  experimental  data  and 
simulations  with  different 
bending  is  shown  in  Figure  6. 

The  figure  shows  that  the 
simulations  with  bending  alone 
do  not  fit  the  experimental  data 
well,  which  indicates  that 
bending  is  not  the  only  factor 
which  causes  the  high  energy 
transmission  to  drop.  However, 
from  Figure  6,  we  can  see  that 
the  range  of  the  bending  radius 
must  be  larger  than  100m  to  give 
the  observed  transmission  at  the 
highest  energy(80KeV). 

Capillary  surface 
oscillations  with  wavelengths 
shorter  than  the  capillary  length  and  longer  than  the  wavelength  of  the  roughness  are 
called  waviness.  The  detailed  shape  of  waviness  is  unknown.  Its  average  effect  can  be 
considered  as  a  random  tilt  of  the  glass  wall,  so  that  the  grazing  angle  of  the  photon  is 
changed  by  a  random  amount,  86,  after  every  bounce.  56  is  a  random  number  between 
-A6max  and  A6max  if  6  >A6max.  The  maximum  random  tilt  angle  A6max  is  an  adjustable 
parameter  which  depends  on  the  waviness  of  the  polycapillary  fiber  To  keep  6'  positive, 
J0is  taken  to  be  a  random  number  between  -6  and  A6max  when  6<A6max .  Since  a 
photon  with  a  incident  angle  smaller  than  A6max  has  a  larger  chance  to  experience  an 

angle  increase  than  an  angle 
decrease,  this  is  physically 
reasonable.  In  Figure  8, 
simulations  with  waviness 
corrections  with  A0max  set  at  1 
mrad  and  2  mrad,  which  is 
comparable  to  the  critical 
angle,  are  compared  with  the 
experimental  data.  This  figure 
shows  that  simulations  with 
Figure  11.  Scheme  of  an  x-ray  beam  interacts  with  a  waviness  alone  do  not  fit  the 
random  tilt  surface.  experimental  data.  This  is 

because  the  waviness 
correction  changes  the 


Figure  10.  Simulated  transmission  spectra  with  different 
bending  and  waviness  compared  with  the  experimental  data  in 
search  for  the  best  fitting  of  Fiber  C. 
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Figure  12.  Three  surfaces,  OA1,  OA2  and  OA3,  with 
different  tilt  angles,  </>l,  <j>2  and  <?>3,  respectively,  from  the 
nominal  surface  OO 


reflected  angle,  not  the  profile.  In  fact 
the  capillary  is  still  considered  to  be 
straight,  so  those  photons  which  have 
few  reflections  will  not  be  significantly 
effected  by  waviness. 

Finally  the  waviness  and 
bending  are  combined  by  increasing 
the  bending  radius  R,  roughly 
determined  in  Figure  6,  and  adding  a 
waviness  parameter,  AGmax.  Several 
trials  are  shown  in  Figure  10.  Sim2  has 
too  much  bending  and  not  enough 
waviness;  sim3  has  too  much  waviness 
and  not  enough  bending;  siml  is  the 
best  fit.  Roughness  is  also  included  in 
those  simulations.  The  source  scan 
simulation  with  the  three  fixed 


parameters  are  plotted  along  with  the  experimental  data  in  Figure  9  for  four  more 
photon  energies.  They  all  fit  quite  well.  This  model,  hereafter  labeled  Ml,  assumed  a 
uniform  distribution  of  tilt  angles,  which  is  not  quite  physical.  In  more  recent  work, 
discussed  as  Model  M2,  it  is  assumed  that  these  tilt  angles,  <{>,  are  normally  distributed 
in  the  range  (-tc/2 ,  n!2)  with  the  mean  value  equal  to  zero.  For  high  quality  optics,  the 
standard  deviation  of  this  normal  distribution,  a,  is  much  smaller  than  the  critical  angle, 


0C.  The  probability  distribution  of  tilt  angles,  <j>,  is 


1 


(1) 


In  this  work,  Model  M2,  consideration  was  taken  of  the  fact  that  the  surface  tilt  angle  will 
affect  the  probability  of  x-ray  impact  on  that  surface.  Taken  to  extremes,  a  surface 
region  perpendicular  to  the  beam  is  much  more  likely  to  intercept  the  beam  than  a 
surface  region  parallel  to  the  beam.  Figure  12  displays  three  surfaces,  OA^  OA2  and 
OA3,  with  different  tilt  angles,  <h,  <t>2  and  <t>3,  respectively,  from  the  nominal  surface  OO’. 
The  projections  onto  the  nominal  surface  for  the  three  surfaces  are  equal,  OA^cos^  = 
OA2*cos<t>2  =  OA3*cos<|>3=  OO’.  The  probability  of  incidence  of  a  parallel  x-ray  beam  with 
incident  angle  0  (with  respect  to  the  nominal  surface  OO’)  hitting  these  tilted  surfaces  is 
given  by  their  corresponding  perpendicular  length  AjBj  (for  j  =  1 , 2  or  3) 

P.  ccAiBl  =  OA,  ■  sin((9  +  </>,)  =  ^  -  ■  sin(<9  +  <j>, ),  ^ 

J  COS  </)j 

We  call  this  the  tilt-corrected  probability  distribution.  The  complete  description  of  this 
distribution  is 
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#(*)  = 


F  sin(  $  +  ,/,) 
cos(  (j,) 

0, 


~0<*K T 


(3) 


where  0  is  the  incident  angle,  4>  is  the  tilt  angle  and  F  is  a  normalization  constant. 

Combining  the  normal  distribution  G(<|>)  with  the  tilt  correction  H(<|>)  gives  <)>  for  a 
certain  incident  angle  0  as  follows: 


m= 


KP  sin(0+^ 
cos^ 


where  0  is  the  incident  angle,  §  is  the  tilt  angle  and  K 


2  , 

—  <(f)<-0 
2 

is  a  normalization  constant. 


Noting  that 

COS (j) 


(5) 


and  that  both  the  incident  angle,  0,  and  the  tilt  angle,  <J>,  are  very  much  less  than  40 
mrad,  we  use  the  approximation 


zt i 

Ce2*1  *(0 +  </>), 


-6  <</)<- 
2  , 

2 


where  the  normalization  constant ,  C,  is 


spin " f?cr  +  2 £J  2 

to  speed  the  calculation.  P(<]))  is  called  the  tilt-corrected  normal  distribution.  In  model 
M2,  a  is  the  variable  that  describes  the  amount  of  waviness.  An  example  of  J(<|>)  and 
P(<t>)  is  shown  in  Figure  13.  The  result  of  the  waviness  calculation  is  shown  in  Figure  14. 
The  model  M2  shows  a  much  larger  drop  at  high  photon  energies  for  lower  values  of 
the  waviness  than  does  model  Ml . 
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angle  <f  when  incident  angle  9  - 
0.009rad,  standard  deviation  a  = 
0.0004rad.  </>  is  in  rad. 
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Figure  14.  Simulations  of 
transmission  spectra  for  fiber  C  with 
only  waviness  compared  with  the 
experimental  data.  This  figure  shows 
the  effects  of  waviness.  It  also  shows 
simulations  using  waviness  alone  do 
not  fit  the  data.  The  simulations  do 
not  include  roughness  or  bending. 


A  comparison  of  the  simulation  with  bending  and 
waviness  to  the  experimental  data  is  shown  in 
Figure  15. 

Using  a  similar  technique,  the  simulation 
was  compared  to  the  fibers  listed  in  Table  1, 
using  the  parameters  listed  in  Table  2.  The 
results  are  plotted  in  Figure  16  and  Figure  17.  In 
Figure  16,  the  transmission  for  fiber  A  shows  a 
rapid  drop  for  energies  above  30  keV.  Although 
fiber  D  has  lower  fractional  open  area  than  fiber 
A,  its  transmission  exceeds  that  of  fiber  A  at 
energies  above  30  keV.  This  is  because  fiber  A 
is  thin  (0.5  mm  in  outer  diameter)  and  flexible, 
therefore  difficult  to  keep  straight  in  the 
measurement  apparatus.  Model  M2,  which  can 
vary  only  bending  and  waviness,  requires  a  much 
sharper  bend  for  fiber  A  than  for  fiber  D.  This 
bending  is  more  significant  at  high  energies, 
where  the  critical  angles  are  smaller. 


Figure  15.  Simulated  transmission  spectra  using 
model  M2  with  fixed  waviness  (<j=0.2  mrad)  and 
different  bending  (bisection  method),  compared 
with  the  experimental  data.  The  simulation  curves 
shown  are  just  a  few  representative  cases  in  the 
bisection  process.  For  M2,  the  best-fit  parameters 
are  cr=  0.2  mrad  and  R  =  225  m;  for  Ml,  they  are  z 
-0.5  nm,  s  =  6  pm,  A6max  =  0.35  mrad  and  R  = 
125  m. 
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Transmission  curves  for  fiber  C  and  D 
are  similar  in  shape  and  are  nearly  flat 
up  to  60  keV,  as  shown  in  Figure  17. 
However,  the  smaller  channel  size  for 
fiber  D  results  in  the  simulation  being 
less  sensitive  to  bending.  Therefore  a 
larger  bending  curvature  (smaller 
radius)  is  required  for  fiber  D.  The 
waviness  correction  for  fiber  D  must 
therefore  be  smaller  to  keep  the 
transmission  about  the  same  as  for 
fiber  C. 


Photon  energy  (KeV) 

Figure  16.  Simulations  of  transmission  spectra 
of  fiber  A  and  fiber  D  with  their  best-fit 
parameters  compared  with  the  experimental 
data.  All  parameters  related  to  this  figure  are 
listed  in  Table  1. 


Fiber  E  is  the  thinnest  fiber. 
Table  2  shows  that  it  has  the  largest 
bending  correction  because  of  its 
flexibility.  Its  transmission  curve  is  flat 
up  to  40  keV.  This  is  due  to  its  very 
small  channel  size,  as  discussed  in  the 
comparison  of  fiber  C  and  fiber  D. 


However,  if  the  channel  size  is  too  small,  as  in  fiber  E,  it  also  results  in  more  reflections 


being  needed  for  a  photon  to  traverse  the  fiber  and  may  have  introduced  other  defects 


such  as  blocked  channels.  This  is  why  the  transmission  is  only  40%  for  the  energies 


below  40  keV  although  the  open  area  is  around  55%. 


Fiber  C:  X  data,-  -  Ml :  best  fit*  A  M2:cM).2,R=225 


Photon  energy  (KeV) 

Figure  17.  Simulations  of  transmission  spectra  of 
fiber  C,  fiber  D  and  fiber  E  with  their  best-fit 
parameters  compared  with  the  experimental  data. 
All  parameters  related  to  this  figure  are  listed  in 
Table  2. 


The  new  model  yields  as  high  a 
quality  fit  with  two  fewer  fitting 
parameters  than  used  in  model 
Ml  in  Figure  4.  The  high 
transmission  and  the  simulation 
results  show  that  the  quality  of  the 
capillary  fibers  is  quite  good.  The 
bending  radius  is  above  130 
meters  for  type  3  and  type  4 
capillaries.  It  is  hopeful  that  we 
can  further  improve  the  high- 
energy  transmission  performance 
of  polycapillary  fibers  by 
decreasing  to  an  optimum 
channel  size,  and  making  them 
more  rigid. 


Ola#1 
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Fiber  Description 

Model  Ml 

Model  M2 

Fiber 

# 

Type 

Outer 

Diameter 

(mm) 

Channel 

Size 

(pm) 

s 

|iim 

R 

m 

A0max 

mrad 

i 

<7 

Mrad 

A 

1 

0.5 

12 

6 

105 

0.4 

1551 

WSEEm 

C 

3 

0.75 

22 

50% 

6 

125 

0.35 

225 

0.2 

D 

4 

4 

12 

55% 

Ha 

6 

110 

0.285 

139 

0.125 

E 

5 

0.3 

K5J| 

55% 

105 

0.7 

6 

28 

0.2 

31 

0.09 

F 

4 

4 

12 

55% 

130 

0.8 

6 

90 

0.45 

90 

0.18 

Table  2  Parameters  for  best-fit  simulations.  R  is  the  bending  radius.  For  Ml,  A0maxis  the 
amount  of  waviness,  z  is  the  roughness  height  and  s  is  the  roughness  correlation  length. 
For  M2,  cris  the  standard  deviation  of  the  waviness. 


Fiber  B  in  Figure  4  has  poor  transmission;  it  is  a  lead  glass  fiber,  with  high 
surface  roughness  and  waviness.  More  recently,  higher  quality  lead  glass  fibers  have 
been  produced.  The  transmission  as  a  function  of  photon  energy  for  these  fibers  is 
shown  in  Figure  18.  These  leaded  fibers  have  higher  absorption  of  high  angle  photons 
than  the  lower  density  borosilicate  glass  fibers,  and  therefore  would  provide  even  better 
scatter  rejection,  as  is  described  in  section  6.3.3. 


Transmission  was  also 
measured  of  deliberately  curved  fibers, 
which  can  be  used  to  deflect  x-ray 
beams  over  larger  angles.  The  results 
are  shown  in  Figure  19. 


x  Experiment 
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a  T  io  T  12  '  14  T  ie  Figure  19.  Transmission  efficiency  of 350  pm  OD 

polycapillary  fiber  with  6  pm  channel  diameter.  The 
fiber  length  is  25  cm.  Simulated  values  are  shown  for 
20,  30,  40,  50  keVxrays.  The  parameters  are 
curvature,  R;  waviness,  w;  and  roughness,  z. 


Figure  18.  Recently  manufactured  lead  glass 
fibers. 


counts/s 
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aperture  location  FWHM(mrad)  peak  center, mrad  6.2.  SIMULATION  OF 

MULTIFIBER  OPTICS 

Measurements  of 
the  exit  divergence  of  the 
collimating  optic  have 
been  performed  at  8  keV 
and  at  20  keV.  Exit 
divergence  is  important  if 
a  capillary  optic  is  to  be 
used  as  a  collector  fore 
slit  before  the  patient.  In 

-0.4  -0.2  0.0  0.2  0.4  xu  I 

this  case,  the  spatial 
resolution  will  depend  on 

Figure  20.  Measured  local  divergence  of  the  output  of  multi-fiber  I  at  8  angular  spread  of 

keV-  x-rays  at  the  exit  end  of 

the  polycapillary  fibers. 

Large  angular  divergence  is  equivalent  to  a  large  focal  spot  size  in  a  conventional 
system. 

The  divergence  can  also  be  used  to  assess  fiber  alignment  in  the  lens.  The 
beam  exiting  the  optic  is  characterized  by  a  divergence  that  arises  from  two  effects.  The 
first  is  the  local  divergence  of  the  x-rays  emerging  from  each  channel,  usually  between 
0C,  the  critical  angle  for  total  reflection  and  20c  There  is  also  fiber  misalignment,  the 
deviation  of  the  individual  channel  axis  direction  from  the  optic  axis  direction.  Local 
divergence  depends  upon  the  surface  roughness,13  waviness,  and  bending  of  fibers 
through  the  optic.  Figure  20  shows  the  measured  local  divergence  of  the  output  of  the 
collimating  lens  at  8  keVwith  the  aperture  at -15  cm,  -10  cm,  0  cm,  10  cm  and  15  cm 
off  the  axis  of  the  lens.  The  FWHM  of  each  divergence  curve  is  in  the  range  of  3.8  mrad 
to  3.9  mrad,  which  is  very  close  to  the  critical  angle  of  about  4  mrad  at  8  kev.  The  peak 
centers  of  these  five  local  divergence  curves  are  systematically  shifted.  This  shift  could 
be  caused  by  the  output  ends  of  the  fibers  inside  the  lens  being  slightly  convergent 

rather  than  parallel.  The  widths 
of  the  measured  local  divergence 
curves  shown  in  Figure  20  are 
larger  than  that  predicted  by  an 
ideal  lens  divergence  simulation, 
shown  in  Figure  21. 

The  divergence  of  the  modeled 
ideal  lens  is  low  because  the 
nearly  straight  central  fibers,  if 
ideal,  would  not  increase  the 
divergence  above  the  entrance 
Divergence  angle  (mrad)  divergence  due  to  the  source 

Figure  21.  An  ideal  whole  lens  simulation  at  8  keV gives  a  Spot  Size.  Figure  23  ShOWS  the 
divergence  FWHM  of  2.5  mrad.  simulated  divergence  profile  of  x- 
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Angle  (mrad) 

Figure  23.  For  an  ideal  fiber  (with  no  surface 
defects),  the  divergence  simulation  yields  a 
FWHM  of  2.36  mrad  at  8  keV. 


Angle  (mrad) 

Figure  22.  A  simulation  with  a  waviness  of  0.15 
mrad  increases  the  divergence  value  from  2.4 
mrad for  an  ideal  straight  fiber  to  4. 0  mrad. 

rays  exiting  from  an  ideal  straight  fiber, 
which  has  a  FWHM  of  2.36  mrad. 


The  simulations  in  Figure  21  and  Figure  23  did  not  include  the  effects  of  profile 
defects  and  waviness.  Waviness  will  increase  the  angle  of  reflection  for  x-ray  photons 
for  most  bounces  inside  the  channel.  Consequently  the  divergence  from  the  lens 
increases.  Figure  22  shows  a  simulated  transmission  of  a  straight  fiber  as  a  function  of 
exiting  angle,  using  a  waviness  of  ©  =  0.15  mrad.  The  width  of  the  simulated 
divergence  curve  for  this  straight  fiber  is  about  4.0  mrad,  which  is  close  to  the  local 


divergence  at  the  center  of  the  lens  shown 
in  figure  16.  Waviness  of  the  channels 
changes  the  divergence  of  the  x  rays 
exiting  from  the  center  of  the  lens  from 
about  2.4  mrad  to  4.0  mrad. 

6.3.  MAGNIFYING  LINEAR 
MONOLITHIC  ANTISCATTER  OPTICS 

6.3.1.  Transmission 

Linear  monolithic  optics  taper  as  a  single 
glass  unit  from  input  to  output  as  shown  in 
Figure  28.  These  optics  are  relatively 
difficult  to  manufacture.  Many  early 
prototypes  have  poor  transmission,  as 
shown  in  Figure  24.  To  help  understand 
the  defect  problem  an  extensive  simulation 
and  measurement  study  was  performed  on 
several  of  the  prototypes. 


Figure  24.  2-D  source  scan  of  early  prototype 
optic  (input  diameter  =  2  mm,  output  diameter 
=  3  mm,  transmission  =  5  %>  at  20  keV,  focal 
distance  22  cm)  This  lens  transmits  poorly 


because  only  the  central  capillary  channels 
Transmission  of  a  monolithic  taper  is  transmit  while  the  outer  capillaries  do  not 
plotted  versus  source  lens  distance  in  transmit  at  all.  Only  about  33  %  of  the  optic 

Figure  25  along  with  the  simulation  values,  cross  section  transmits  at  5  %  or  greater. 
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250  450  650  850  Photon  energy  (keV) 


Source  lens  distance  (mm)  Figure  26.  Transmission  of  whole  tapered  lens. 


Figure  25.  Transmission  of  tapered  lens. 
Solid  line  is  simulation. 


Pinhole  size:  1 


Pinhole  size:  0.5 


Photon  energy  (keV) 

Figure  27.  Transmission  of  the  central  part  of  the 
tapered  optic. 


The  maximum  point  for  transmission  is  j=  30  -  *  ♦  mm 

around  450  mm.  The  second  method  s  20  -  — simulated 

can  be  affected  by  the  source  "  10-  i“w^ness 

inStabi'ity-  .  .  0  20  40  60  80  id 

Transmission  was  measured  in  photon  energy  (keV) 

a  wide  energy  range  for  the  whole  lens 

as  well  as  the  central  part  Of  the  lens.  FiSure  27 ■  Transmission  of  the  central  part  of  the 

The  transmission  spectrum  for  the  tapered  optic. 

whole  lens  is  plotted  in  Figure  26.  The 

figure  shows  that  the  transmission  drops  quickly  for  photons  with  energy  higher  than  20 
keV,  but  the  transmission  under  20  keV  is  almost  flat.  The  low  transmission  in  the  high 
energy  range  can  be  explained  by  the  non-linearity  of  the  lens.  The  outermost  channels 
are  more  bent  than  the  central  channels,  so  they  have  low  or  zero  transmission  at  the 
higher  photon  energies.  This  is  confirmed  by  the  higher  transmission  measurement  at 
the  central  part  of  the  lens,  shown  in  Figure  27  and  the  transmission  uniformity 

measurement  at 

aft  different  energies, 

fore  phantom  e0iiimatio_n  shown  in  Figure 

c°",nnati^n  ■  29.  Transmission 

Itj^m  k  k  ||  uniformity  of  the 

. - z=Mz= |B  lens  was  measured 

if  I  1 1  by  putting  a  200 

source  T  II  pm  pinhole  in  front 

optic  ®  of  the  lens  and 

detect/  performing  a  two 

dimensional  scan. 

Figure  28.  Small  optic  scanning  geometry  for  contrast  and  imaging  At  8  keV,  the  whole 

measurements.  |ens  transmits, 


fore 

collimatipn 


aft 

collimation 


Figure  28.  Small  optic  scanning  geometry  for  contrast  and  imaging 


measurements. 
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Figure  29.  Two  dimensional  scans  of  the  output  of  the  tapered  optic. 

although  we  can  see  the  transmission  is  not  completely  flat  due  to  the  non-linearity  and 
the  defects  of  the  lens,  such  as  blocked  channels.  At  25  keV,  transmission  of  the  whole 
lens  drops  to  25%.  Most  of  the  lens  is  still  transmitting  at  25  keV,  but  the  outer 
channels  have  lower  transmission  than  the  central  part.  At  35  keV,  the  lens  has  a 
narrower  transmitting  region.  The  measurements  show  that  the  lens  has  undesirable 
nonlinearity. 

6.3.2.  Simulations 

Modeling  defects,  including  waviness,  bending,  and  blockages,  require  sophisticated 
simulations.  The  simulations  plotted  in  Figure  25  and  Figure  26  used  the  same 
roughness  and  waviness  correction  as  that  used  for  fiber  D  in  Figure  4.  A  linear 
tapered  lens  is  ideal  for  mammography,  because,  as  shown  above,  any  curvature  can 
cause  transmission  reductions.  However,  the  lenses  currently  available  have  nonlinear 
profiles.  The  profile  is  not  circular,  but  to  simplify  the  simulation,  we  assume  that  the 
profile  is  circular.  A  uniform  bending  will  provide  the  best  transmission  performance  for 


Energy  (keV) 

Figure  30.  Transmission  of final  half  of  A  (A2),  and 
its  two  constituent  pieces ,  the  first  5.5  cm,  A21  and 
the  final  5.5  cm,  A22.  A21  contains  a  defect. 


Energy  (keV) 


Figure  31.  Transmission  of  original  optic  compared 
to  two  pieces.  B  is  the  first  10  cm  (closest  to  the 
source),  A  is  the  final  22  cm. 
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a  given  bending  angle,  so 
the  simulation  will 
overestimate  the  lens 
performance.  A  further 
simplification  assumes  the 
lens  has  a  fixed  channel 
size,  although  the  channels 
actually  taper  with  the  lens. 

The  geometry  for 
uniformly  bent  nonlinear 
lenses  is  shown  in  Figure 
32.  The  simulation  method 
for  a  uniformly  bent  lens  is 
an  extension  of  the 


Figure  33.  Curvature  estimation  for  lens  simulation. 


simulation  for 

uniformly 

bent 

polycapillary 
fibers.  The 
simulation  is 
done  by 
sampling 
channels  in 


the  input 

Cross  section  figure  32.  Geometry  for  Lens  simulation. 


of  the  lens 


with  an  even  step  size.  In  the  simulation  for  a  bent  fiber  it  is  assumed  that  the  capillary 
bends  in  -y  direction,  and  z  direction  is  along  the  input  axis  of  the  fiber  In  order  to  reuse 
the  old  simulation  code,  the  coordinate  of  the  simulation  x-y-z  is  rotated  to  x’-y’-z’, 
where  z’  is  the  direction  of  channel  axis  at  the  input  end  and  -y’  is  the  channel  bending 
direction.  Assuming  the  source  position  and  the  capillary  position  at  input  end  is  (xs\  ys’, 
zs’)  and  (0,  0,  Zj’)  respectively,  then  the  source  capillary  distance  is  z--  zs’.  The  origin  of 
the  coordinate  is  selected  to  be  the  focus  point  of  the  lens. 

Another  parameter  needed  by  the  simulation  is  the  bending  curvature  of  the 
selected  channel.  This  value  varies  with  the  position  of  the  channel.  The  calculation  of 
the  bending  curvature  is  demonstrated  in  Figure  33. 

First  it  is  assumed  that  the  channels  are  evenly  distributed  on  both  input  and 
output  cross  section.  Let  input  and  output  diameter  of  the  lens  be  r and  r0,  and  focal 
distance  ~oa  be  f,  and  the  lens  length  ac  be  L,  then 


AB 


For  a  =  AB  /  CD  d  =  CE  -  CD 

a  «  L  /  R  =  d  /  (L  /  2) 
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where  R  is  the  bending  radius,  d  is  the  distance  between  point  D  and  E  as  shown  in 
Figure  33,  and  a  is  the  bending  angle.  Thus 

R  =  L2  /  2  d  (8) 

The  simulated  results  in  Figure  26  and  Figure  27  show  higher  transmission  than 
the  experimental  results.  This  indicates  that  the  lens  has  more  reduction  in  transmission 
due  to  the  bending.  This  is  expected,  because  the  lens  does  not  have  a  uniform 
bending. 

6.3.3.  Potential  for  Lead  Glass 

Transmission  from  early  fibers  was  poor  in  the  whole  energy  range.  Recent 
fibers  are  more  promising.  Transmission  for  an  early  95  mm  long  and  a  more  recent 
125  mm  long  lead  glass  polycapillary  fiber  is  shown  in  Figure  34.  The  transmission 
drop  at  lower  energies  may  be  caused  by  partially  blocked  channels.  Transmission  can 
also  be  improved  by  using  a  short  optic.  Almost  total  scatter  rejection  is  expected  for  a 
lead  glass  optic  longer  than  30  mm. 

6.3.3. 1 .  Calculation  of  Contrast  Improvement  and  Relative  SNR 
Lenses  with  large  diameters  are  not  currently  available.  Their  potential  behavior  can  be 
predicted  by  scaling  up  existing  measurements  to  a  large  area  lens,  as  shown  in  Figure 
35.  The  value  of  h  is  the  length  of  the  optic.  For  a  point  P  on  the  film,  scattered  radiation 
could  come  from  any  direction  in  the  solid  angle  Qpatien„  which  is  the  solid  angle 
subtended  by  the  patient.  The  transmission  for  scattered  x  rays  can  be  calculated  from 
the  mass  aborption  coefficients.  The  calculated  scatter  transmission  depends  on  the 
length  of  the  photon  path  through  the  optic.  Making  a  very  conservative  estimation,  we 
can  fix  this  length  to  the  length  of  the  optic.  There  is  also  a  small  part  of  the  scattered  x 
rays  that  have  almost  the  same  direction  as  those  primary  x  rays.  This  part  of  the 
scattered  x  rays  with  an  incident  angle  less  than  critical  angle  could  still  pass  through 
the  channel.  So  the  total  scatter  transmission  of  the  lens  is 
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Figure  34.  Transmission  as  a 
function  of  photon  energy  fora 
lead  glass  capillary.  The  fibers 
are  105-120  mm  in  length. 


Figure  35.  An  ideal  large  area  lens. 
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Table  III.  Transmissions  for  three  lenses.  The  results  for  the  ideal 
borosilicate  lens  are  measured  results  for  the  center  part  of  the 
prototype  lens.  The  results  for  lead  glass  lens  are  simulated  with 
120  kVp  spectrum  and  a  W  target. _ 


Lens 

# 

Lens  type 

Length 

(mm) 

Photon 

Energy 

Primary 

transmission 

TD-0DtiC 

Scatter 

transmission 

Ts-ODtiC 

1 

Ideal 

borosilicate 
Lens.  Scaled 
up  from  the 
center  part  of 
the  prototype 
lens 

166 

20  KeV 

0.6 

0.0 

45  KeV 

0.55 

0.002 

60  KeV 

0.55 

0.005 

70  KeV 

0.45 

0.03 

2 

Actual 

Borosilicate 

prototype 

lens 

166 

20  KeV 

0.33 

0.0 

3 

Lead  glass 

30 

69  KeV 

0.63 

0.0076 

_  Ns_„plk  _  T  Tp  •  7tdc 

*  s-optic  -  -  1  c 


N. 


Q 


patient 


(9) 


where  Ns.optic  is  the  number  of  scattered  x  rays  with  the  optic,  Ns  \s  the  number  of 
scattered  x  rays  without  the  optic,  Tsis  the  transmission  for  large  angle  scattered  x  rays, 
Tp  is  the  transmission  for  primary  x  rays,  Qpatien,is  the  angle  subtended  by  the  patient,  0C 
is  the  critical  angle  for  transmission  in  the  channels.  Because  the  critical  angle  is  on  the 
order  of  10'3  radians,  the  second  term  is  very  small,  and  can  be  ignored.  So  scatter 
transmission  of  the  lens  can  be  estimated  to  be  equal  to  large  angle  scatter 
transmission,  Ts,  which  can  be  calculated  from  equation  Error!  Reference  source  not 
found.,  and  nearly  scatter  free  lenses  can  be  obtained  by  using  a  suitable  length  optic. 
In  the  above  equation,  we  assumed  that  the  angular  distribution  of  the  scattered  x  rays 
is  isotropic.  If  the  scatter  distribution  is  not  isotropic,  the  second  term  in  equation  (9)  can 
be  larger  than  expected,  but  will  still  be  very  small,  and  this  assumption  will  not 
significantly  effect  the  result.  In  equation  (9),  Ts ,  Tp  and  9C  are  all  energy  dependent.  In 
the  case  where  x  rays  with  a  wide  x-ray  spectrum  are  used,  Ts  and  Tp  should  be 
replaced  by  the  average  large  angle  transmission  and  primary  transmission.  To  avoid 
underestimating  scatter  transmission,  6C  can  be  set  to  be  the  critical  angle  for  the  lowest 
energy  in  the  spectrum. 

The  scatter  and  primary  transmission  of  the  optic  determine  the  performance  of  the 
capillary  lens.  For  example,  if  the  scatter  fraction  before  the  capillary  lens  is  F,  then  the 
scatter  fraction  after  the  lens,  Foptic,  is 
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F -T  ■ 

p  _  s-optic 

op,!c~  F-Ts_oplic+(\-F)-Tp 


1  -F(l-K) 


where  Tp.optiCis  the  primary  transmission  of  the  lens  and  K  is  the  ratio 


The  contrast  enhancement  achieved  by  using  the  optics  is 


FK  +  (l-F ) 


The  contrast  enhancement  alone  is  not  the  only  relevant  quality  factor.  The  signal  to 
noise  ratio  (SNR)  is  another  important  quality  factor.  According  to  the  Rose  model, 
k2  =  C2<&A  ,  (13) 

where  k  is  the  SNR,  C  is  the  contrast,  ®  is  the  photon  flux,  and  A  is  the  area  of  the 
target.  In  order  to  compare  the  SNR  for  different  scatter  rejection  devices,  relative  SNR 
was  used.  Relative  SNR  is  the  SNR  normalized  by  the  SNR  without  any  scatter 
rejection  device, 

*/=7T  =  ^r(FTs+(l-F)Tp),  (14) 

fc0  C0 

where  kr  is  the  relative  SNR.  C  and  C0  are  the  contrast  with  and  without  out  the  scatter 
rejection  device,  respectively.  Ts  is  the  scatter  transmission  and  Tp  is  the  primary 
transmission  of  the  device. 


6. 3. 3.2.  Performance  of  Three  Hypothetical  Lenses 

Performance  estimates  of  three  hypothetical  lenses  are  given  in  Table  III.  The  ideal 
borosilicate  lens  is  scaled  up  from  the  center  part  of  the  measured  lens.  Simulated 
transmission  was  used  for  the  lead  glass  lens.  The  performance  of  the  three  lenses  was 
compared  with  the  performance  of  the  conventional  grids  described  in  Table  IV.  The 
primary  and  scattered  transmissions  of  the  different  grids  in  Table  IV  were  measured  at 
69  keV,  which  is  the  average  energy  of  the  120  kVp  tungsten  spectrum,  and  with 

scatter  fraction  of  85%,  by  C.  E.  Dick  and  J. 


Table  IV.  Performance  for  several 
commercial  grids,  meassured  at  a 
scatter  fraction  of  85%. 14 


Grid  ratio 

Primary 
transmission 
at  69  keV 

Scatter 
transmission 
at  69  keV 

6 

0.73 

0.27 

8 

0.64 

0.13 

12 

0.61 

0.071 

W.  Motz  14.  It  is  assumed  that  the  scatter  and 
primary  transmission  are  the  same  under 
different  scatter  fraction.  This  assumption  is 
correct  if  the  angular  distribution  of  the 
scattered  x  rays  does  not  change.  The 
performance  of  the  grids  is  assumed  to  be 
energy  independent.  They  might  perform 
better  at  lower  photon  energies  if  the  scatter 
transmission  is  due  to  the  lead  plates  which 
are  not  thick  enough  to  absorb  the  incident  x 
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Figure  36.  Contrast  improvement  vs.  scatter 
fraction  for  the  three  lenses.  The  results  were 
compared  with  the  performance  of  two 
commercial  grids. 


rays.  However,  it  was  assumed  the 

10  _ 

O  Simulated  Borosilicate,  20  KeV  dominant  mechanism  for  scatter 

|  9  a  simulated  Borosilicate, 60  Kev  •  transmission  was  x  rays  passing 

ra  8  -  O  Simulated  Borosilicate.70  KeV  .  . 

I  7_  x  Whole  Measured  Lens,  20  Kev  through  the  interstitial  material.  The 

£  •  Simulated  Lead  Glass  a  contrast  enhancement  and  relative 

I  ,  Grid  ratio  6  SNR  were  calculated  under  different 

o_  5  aBaaHBGrid  ratio  1 2  jg) 

|  4  _  scatter  fractions  and  plotted  in 

|  3  Figure  36  and  Figure  37. 

<S  2  -  For  most  medical  x-ray  imaging 

i  1  a-,*  — — I  applications,  the  effective  energy 

o.2  o.3  o.4  o.5  o.6  o.7  o.8  o.9  ranges  from  20  to  70  keV. 

scatter  fraction  Mammography  has  an  effective 

energy  of  20  keV.  The  results  in 

Figure  36.  Contrast  improvement  vs.  scatter  Fi9ure  36  show  that  contrast 

fraction  for  the  three  lenses.  The  results  were  improvement  of  all  the  three 

compared  with  the  performance  of  two  hypothetical  lenses  are  ettert  an 

commercial  grids.  f r  as  9°od  as  hat  of  he  "T3  9"d,s- 

Lens  #2  is  scaled  up  from  the  whole 

measured  lens.  This  lens  has  good  contrast  improvement  because  of  its  total  scatter 
rejection.  However  its  relative  SNR  is  lower  than  that  of  the  grids  when  the  scatter 
fraction  is  less  than  80%  because  of  its  low  primary  transmission.  The  primary 
transmission  of  this  lens  is  only  33%  and  the  primary  transmissions  for  all  the  three 
commercial  grids  are  better  than  60%.  So  with  this  lens  more  exposure  will  be 
necessary  to  the  patient.  Overall  primary  transmission  has  to  be  improved  for  practical 
use.  Lens#1  is  scaled  up  from  the  center  part  of  the  measured  lens.  From  the 

2.5  i -  , ,  — — ,  „„„„ - 3—1  calculations,  summarized  in  Figure  36 

o  Simulated  Borosilicate,20  KeV  %  r 

2.3  A  Simulated  Borosilicate, 6o  Kev -  and  Figure  37,  its  performance  is 

21  o  simulated  Borosiiicte, 70  Kev  _  significantly  better  than  those  of  the 

1*„_  ;  “sLeaf2oGSv'ens  _  commercial  grids  for  effective  energies 

2  ^  y  -  -  «Gridratio6  4  &  I  of  20-45  keV.  This  lens  will  not  only 

”  _  " - Grid  ratio  12  Y  improve  the  contrast  in  the  imaging  but 

~  also  reduce  the  exposure  time  to  the 

^  1.3  "  patient.  So  it  is  ideal  for  mammography. 

Q£  i.i - -  Assuming  that  the  lead  glass 

0.9  __9 -  capillaries  have  the  same  quality  as  that 

0.7  -  -  H  -  x _  of  borosilicate  fibers,  we  can  estimate 

0.5  X  the  performance  of  a  lead  glass  lens. 

02  03  04  05  06  07  08  09  Simulations  of  transmission  spectra 

Scatter  fraction  '  '  were  performed  for  a  source  to  lens 

distance  of  1.5  m,  and  the  source 
operated  at  120  kVp  with  a  0.5  mm  spot 

Figure  37  Relative  Signal  to  Noise(SNR)  vs.  size.  For  a  perfect  lens,  each  channel 

scatter  fraction  for  three  hypothetical  lenses,  in  the  lens  is  tapered.  With  a  tapered 

The  results  were  compared  with  the  channel,  fewer  average  reflections  will 

performance  of  three  commercial  grids.  be  needed  for  the  photons  to  pass.  So 
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Figure  37  Relative  Signal  to  Noise(SNR)  vs. 
scatter  fraction  for  three  hypothetical  lenses. 
The  results  were  compared  with  the 
performance  of  three  commercial  grids. 
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theoretically,  transmission  with  a 
perfect  lens  is  larger  than  for  a 
straight  capillary  if  they  have  the 
same  quality.  The  transmission  is 
quite  sensitive  to  the  capillary  length. 

Therefore,  as  long  as  the  lens  has 
low  enough  scatter  transmission,  the 
length  of  the  lens  should  be  kept  as 
short  as  possible.  The  performance 
estimation  results  of  the  lead  glass 
lens  in  Figure  37  are  based  on  the 
simulated  transmission  with  length  of 
30  mm.  The  scatter  transmission  of 
this  lens  is  less  than  1%.  Primary 
transmission  as  high  as  that  of  grid  is 
also  obtained  with  simulation. 

Compared  to  the  grid  performance,  Figure  39  Three-dimensional  edge  image  after 
this  lens  has  higher  contrast  background  subtraction. 

enhancement  and  relative  SNR.  So 

this  lens  will  not  only  improve  the  contrast  of  the  image,  but  also  reduce  the  radiation 
dose  to  the  patient. 

6.3.4.  MTF  Calculation 

Modulation  Transfer  Function  (MTF)  is  the  most  fundamental  measurement  of 
spatial  resolution  used  in  radiology.  The  standard  technique  is  to  image  a  slit, 
determine  the  line  spread  function  (LSF),  and  compute  the  Fourier  transform.  The  MTF 
measured  for  the  mammography  system  with  capillary  optics,  compared  to  that  without 
the  capillary  optics,  will  give  the  resolution  improvement  and  also  can  be  used  to 
diagnose  artifacts. 

The  structure  of  the  capillary  bundles  was  obvious, 
not  be  obtained  by  imaging  a  slit  without  background 
elimination.  The  background  deducted  slit  image  is  the 
slit  image  divided  by  the  image  with  capillary  only. 

However,  this  involves  image  registration.  The 
information  for  background  that  can  be  used  for  image 
registration  is  too  little  in  the  slit  image.  An  alternative 
method  is  calculating  the  LSF  as  the  derivation  of  the 
edge  spread  function  (ESF).  The  edge  image  left 
enough  background  for  registration  so  that  it  could  be  a 
feasible  method  in  our  case. 

A  registration  algorithm  (FMI-SPOMF)  proposed 
by  Chen,  et  al. 15  was  used  in  this  work.  This  is  a  method 
to  match  a  two-dimensional  image  to  a  translated,  rotated 
and  scaled  image.  The  approach  consists  two  steps:  the 
calculation  of  a  Fourier-Mellin  Invariant  (FMI)  descriptor 


In  this  case,  the  LSF  could 


Figure  38.  Magnified  image  of  a 
lead  edge.  Real  image  size: 
5(mm)  X  5(mm). 


t.  » 
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Figure  41.  Reprojection  of  a  two- 

dimensional  edge  image  into  a  finely  Figure  40  Edge  spread  function  and  its  fitting  curve, 

sampled  ESF. 

for  each  image  to  be  matched,  and 

then  matching  of  the  FMI  descriptors.  FMI  descriptors  are  translation  invariant.  The 
matching  of  the  FMI  descriptors  is  to  find  out  the  rotation  and  scaling,  and  achieved 
using  a  symmetric  phase-only  matched  filtering  (SPOMF).  In  our  case,  there  is  no 
change  in  scaling,  but  rotation  is  probably  involved.  When  the  rotation  was  found, 
image  translation  is  found  by  SPOMF  method. 

Images  were  taken  with  and  without  edge  as  shown  in  Figure  38.  The  edge  was 
made  of  a  lead  plate.  The  size  of  each  image  is  fifty  by  fifty  in  pixels.  They  were  taken 
in  an  experimental  digital  mammography  system  with  a  computed  radiography  digital 
phosphor  plate.  The  digital  CR  plate  has  limited  resolution,  5  Ip/mm,  and  is  usually 
considered  to  be  not  good  enough  for  clinical  mammography.  But  the  effective 
resolution  of  the  system  could  be  improved  by  the  magnification  of  the  capillary  optics. 
The  magnification  of  the  capillary  lens  used  in  the  system  is  1 .86,  so  that  the  effective 
resolution  should  be  improved  by  a  factor  of  1.86. 

The  registration  algorithm  was  implemented  in  IDL.  With  the  FMI-SPOMF  registration 
algorithm,  no  rotation  was  found.  The  translation  was  then  found  by  SPOMF  algorithm. 
The  difference  between  the  usual  correlation  method  and  SPOMF  method  is  that  it  only 
uses  the  phase  information.  The  phase-only  correlation  function  has  sharper  peak  than 
normal  correlation  function.  Sub-pixel  resolution  could  be  achieved  theoretically. 
Flowever  it  was  found  that  the  resolution  is  not  good  enough  in  our  case.  This  may  be 
caused  by  the  relatively  small  image  size.  The  background-subtracted  result  was  further 
optimized  by  manual  shifts  in  sub-pixel  range.  The  result  after  manual  optimization  was 
shown  in  Figure  39  .  The  periodic  background  is  totally  gone.  A  cubic  spline 
interpolation  method  was  used  for  the  best  result  while  shifting  the  background  image 
relative  to  the  edge  image. 
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Figure  42  Calculated  MTF  compared  with  the 
MTFs  with  scan  optics,  and  without  optics. 


The  presampling  MTF  includes 
unsharpness  of  the  detector  and  the 
sampling  aperture.  In  order  to  eliminate  the 
aliasing,  a  finely  sampled  ESF  is  obtained 
with  a  slightly  angulated  edge  in  a  single 
exposure  as  illustrated  in  Figure  41. 16,17  The 
angle  in  our  edge  image  in  Figure  39  was 
calculated  to  be  around  6  °.  The  resultant 
ESF  is  plotted  in  Figure  40  .  A  direct 
calculation  of  MTF  is  to  take  the  derivation 
of  the  ESF  and  result  a  LSF,  then  take 
Fourier  transform.  As  it  can  been  seen,  the 
ESF  in  our  case  is  not  smooth  enough  for  a 
direct  calculation.  An  alternative  technique 
is  to  use  an  ESF  fitting  procedure.18  In  the 
fitting  method,  the  ESF  is  represented  with 
a  term  (1-exp)  and  an  error  function  (erf)  as 
in  equation  (15).  MFT  is  calculated  by 
equation  (16),  where  a,  b,  c  and  d  are  fitting 
parameters  from  equation  (15).  The  fitting 


curve  is  plotted  in  Figure  40  as  the  solid  line. 


£,SF(x)  =  a{l-exp(-&|x-x0  | )}  +  c  •  erf  (dil2  |x-x0 1)  O5) 


MTF(f)  = 


cex  p(-^2/2/^)  +  a(l  +  47r2/2/62)~' 
(c  +  a) 


(16) 


The  calculated  MTF  was  compared  with  that  from  scanned  optics  in  Figure  42  .  The 
MTF  from  scanned  optics  was  calculated  from  a  slit  image.  Optic  structure  was 
smeared  out  with  scanning,  so  no  background  deduction  is  necessary.  An  ideal  MTF 
was  also  plotted  in  Figure  42  .  An  ideal  MTF  is  the  one  when  we  assume  the  lens 


Amptek  detector 


X-Y -Z  stages 


Figure  43.  Schematic  of  the  experimental  setup  for  optic  characterization. 
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introduces  only  the  1.86  factor  magnification  and  no  image  degradation.  It  is  measured 
with  a  slit  image  without  optics,  and  its  result  was  multiplied  by  the  magnification  of  the 
capillary  lens.  The  measured  MTF  with  lens  should  no  better  than  the  ideal  MTF. 
However,  for  some  reason,  the  MTF  for  stationary  lens  is  better  than  the  ideal  MTF.  It  is 
possible  that  the  background  elimination  process  failed  the  MTF  measurement.  Another 
possibility  is  that  the  magnification  of  the  lens  is  different  from  1 .86,  since  the 
magnification  was  measured  separately  at  a  different  time. 
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7.  Key  Research  Accomplishments 

•  Development  of  extensive  modeling  capability  for  polycapillary  optics  that  includes 
realistic  models  for  the  effects  of  profile  defects  and  shows  extremely  good 
agreement  with  measured  data  and  therefore  the  capability  to  predict  performance 
in  new  geometries. 

•  Measured  polycapillary  fiber  transmission  in  excess  of  65%  at  20  keV  for 
borosilicate  glass  and  in  excess  of  50%  for  lead  glass.  Measured  high  angle 
transmission  (which  corresponds  to  scatter  transmission)  is  lower  than  1%  for  12  cm 
long  borosilicate  glass  optics  and  1  cm  long  lead  glass  optics. 

•  Modeling,  measurement  and  defect  studies  have  led  to  better  manufacturing 
capabilities  for  magnifying  tapered  monolithic  optics.  These  optics  have 
demonstrated  the  ability  to  provide  nearly  ideal  contrast  enhancement,  while  at  the 
same  time  increasing  the  system  MTF  performance  at  all  spatial  frequencies.  A 
multitaper  optic  has  been  produced  and  tested. 

•  Image  analysis  studies  have  provided  the  means  for  image  artifact  analysis  and 
reduction. 
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9.  Conclusions 

Some  manufacturing  problems  have  occurred  with  scaling  up  the  optics  to 
clinical  size.  These  problems  have  been  addressed  by  the  studies  of  defect  properties. 
Notwithstanding  the  problems,  the  optics  tested  have  performed  extremely  well. 
Multifiber  collimating  optics  and  monolithic  linear  magnifying  tapers  are  being  studied. 
Scatter  rejection  is  very  high  from  all  the  optics.  The  linear  tapers  resulted  in  nearly 
ideal  contrast  enhancements  and  simultaneously  increase  in  MTF  at  all  spatial 
frequencies.  This  will  result  in  an  improvement  in  resolution  regardless  the  inherent 
spatial  resolution  of  the  detector.  The  rapidly  growing  modeling  capability  has  lead  to  a 
real  development  of  understanding  of  the  nature  of  polycapillary  defects,  which  is 
already  leading  to  improvements  in  the  manufacturing  processes.  Further,  the 
increased  confidence  in  the  modeling  allows  future  lens  geometries  to  be  effectively 
“tested”  in  simulation  so  that  design  parameters  can  be  rapidly  optimized. 
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